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Discovery of the Higgs boson in any decay channel depends on the existence of event variables or 
cuts with sensitivity to the presence of the Higgs. We demonstrate the non-optimality of the kine- 
matic variables which are currently expected to play the largest role in the discovery (or exclusion) 
of the Higgs at the LHC in the rr channel. Any LHC collaboration looking for opportunities to 
gain advantages over its rivals should, perhaps, consider the alternative strategy we propose. 



Introduction 

There is much to be gained from constructing event 
variables which place maximal lower-bounds on well de- 
fined quantities of interest. Such variables can be used 
to select events containing new-physics when the scale of 
the property which is being "bounded" is higher in the 
signal than in the most important backgrounds. One 
may construct the single variable that bounds an arbi- 
trary scale by considering that scale (often a mass) to be 
a function of all the unknowns in the event (often compo- 
nents of invisible particle momenta). Having done this, 
the minimal value of this scale over all possible values of 
those unknowns, subject to any constraints that need to 
be asserted to enforce consistency, is the bound in ques- 
tion.^ The transverse mass is an example of such a max- 
imal lower-bound variable: when applied to a W ^ lu 
event in a hadron collider it returns the largest possible 
lower-bound on the M^-mass that may be derived from 
that event (given access to the lepton four-momentum 
and the missing transverse two-momentum only) assum- 
ing that there were no confounding sources of missing 
transverse momentum. 

Separating Z ^ tt from h ^ tt 

The main background to searches for h ^ tt is 
Z TT, so following the general procedure described 
above, one would expect that the best way to separate 
the signal from this irreducible background is to con- 
struct the variable which provides the maximal lower 
bound for the "parent" mass (i.e. mjj or mz in signal 
and background respectively) given the observed visible 
decay products of the taus together with the net miss- 
ing transverse momentum. In a "perfect" detector, such 
a variable (we will call it TT^H^sgs-bound-j g^ould, by con- 
struction, place all the irreducible background Z — >■ rr 
events at values of rn^^^ss -hound j^g^Q-^^ j^^^ leaving the 



^ Examples in the context of Higgs boson searches include Refs. [J 
[2 . In Ref . Q a recent attempt has been made to formally write 
down the steps that are needed to construct such maximal lower- 
bounding variables for a wide class of circumstances. 
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FIG. 1: At Monte Carlo truth level, and for signal only, the 
distribution of the new m^.^^^""''"""'^ variable (black) is com- 
pared to the existing variables m"^^" [l|] (magenta), m^f'''^*'^'^ 
11 (blue) and my^"''''" [1, 0] (cyan) for a 120 GeV Higgs de- 
caying to TT. All histograms are scaled to unit area. 



region mz < m^^sgs- bound ^ available to the signal 
and entirely free of background. Following the general 
procedure described above, we are therefore naturally led 
to construct TTjH^ggs-bound defined as follows^: 

^Higgs-bound ^ min JWW: (1) 

where 

H^' = PC + Qt + P.^ + (2) 

is the four momentum sum of the measured visible 
and hypothesised invisible 2 momenta of the daugh- 



^ The notation used in this letter follows that of 3] exactly - see 
in particular Tables 1, II, VI and VII therein for reference. In 
the specific context of the decay — >■ T1T2 we denote the mea- 
sured four momentum of the visible decay products of the harder 
and softer tau by and respectively. Each tau has decay 
products (one or more neutrinos) which are unobservable. We 
cannot measure the momenta of these decay products, but we de- 
note hypothesised values for them as and Qj- We have not 
been able to obtain a simple algebraic form for m^|^^'^~''°""'^; 
instead we evaluate it via a computer algorithm. The relevant 
code is available from the authors on request. 
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ters of the two taus, and where and Q2 are subject to 
constraints H comprising: four internal mass constraints 



0, 
0, 



Pf)(QiM + AM) 



(3) 

(4) 
(5) 
(6) 



and one constraint on the missing transverse momentum 
two- vector 



(7) 



Finahy we note that it may be shown that there exists at 
least one pair of momenta Q'^ and Q2 satisfying all the 
constraints if and only if 



Mt2{Pi,P2Jt) < 



(8) 



where Mt2 is the stransverse mass 0-Q ■ Accordingly, it 
is necessary to impose a pre-selection ^ on events before 
^mggs-bound computed. 

Note the difference between the design of TOH^ggs-bound 
and the design of another maximal lower-bound kine- 
matic variable, m!^''"°, ili] which was proposed for Higgs 
mass measurement in h ^ WW — Ivliy events.^ The 
key difference between the tt and the WW topologies 
(other than the obvious fact that m-r <C mw) is that the 
W^s need not be near their mass shells, particularly when 
mil < 2m\Y- Consequently m^^'' does not enforce the 
intermediate W mass-shell constraint.'* 

Simulations 

To compare the performance of Tj^H^ggs-bound against 
other mass-scale variables, we simulate both the signal 
process h ^ tt and the dominant background — tt 
using the HERWIG 6.505 [H, Monte Carlo generator, 
with LHC beam conditions {y/s = 7 TeV). 

The generated tau leptons can decay either leptonically 
(e.g. — > e^Vei^r) or hadronically (e.g. Xv^, 



Note that rriT^^" is exactly the same as the older "cluster trans- 
verse mass" of 0, [Toll which was proposed for the same purpose. 
It is regrettable that a new notation for an existing quantity was 
introduced in [l| , whose authors were sadly not aware of 0, [lol | 
at time of publication. Without prejudice to earlier work, we 
nonetheless retain the ml^^'^ notation to keep a consistent no- 
tation with the papers to which this work is most closely tied: 

Note that one can define a variable that applies only one internal 
W-mass constraint. Such a variable may be better than mrj^''"'^ 
at measuring mj, when mj, < 2myY , however it is unlikely to be 
better for Higgs discovery since there is no resonant background 
of the form Z — > WW — > Ivlu that needs to be suppressed. 
For an example of a variable that is not constructed as a mass- 
bound variable in the context of the h — > WW — lulv channel, 

:i2i. 



where X consists of hadrons or their subsequent decay 
products). The momenta of the visible daughters from 
the tau lepton decays - electrons, muons, hadrons and 
photons - ought to be well-measured by the LHC exper- 
iments. By contrast the contribution of the neutrinos to 
fr^ must be inferred from the negative sum of the mo- 
menta of all observed particles and so can vary consider- 
ably from its ideal value. In our simulations the missing 
transverse momentum is reconstructed as 



E 



PT,i 



where the first sum runs over all reconstructed jets, and 
the second runs over any stable particles within fiducial 
pseudorapidity {\ri\ < 5) and momentum (px > 0.5 GeV) 
that are not clustered into jets. 

The jets used to calculate -^^ are reconstructed using 
the fastjet [ist implementation of the anti-fcx algorithm 
[1^, using the E combination scheme, with distance pa- 
rameter R = 0.6 and minimum jet of 15 GeV. Their 
energies are smeared by a Gaussian probability density 
function of width 



(j{E)/Ej = (o.6GeV5/A 



0.03 



where Ej is the unsmeared jet energy. This resolution is 
ty|)ical of one of the general-purpose LHC detectors [13, 

In this illustrative example, all combinations of 
hadronic and leptonic tau decays are treated on the same 
footing. We select events that contain two taus with 
pseudorapidity satisfying I77I < 2.5, the typical angular 
acceptance of the tracking detector. We require that the 
visible decay products (whether electrons, muons or tau- 
jets) have pr > 20 GeV, and that > 20 GeV. 

As noted earlier, when plotting 777,H^ggs-bound addi- 
tionally require ([8]) to ensure the existence of a minimi- 
sation domain in ([1]). We note that in the narrow-width 
limit, well- measured tau pair events should satisfy dS]) 
by construction. Therefore ui^^sss-bound guarenteed 
to exist in the idealized case. Detector resolution effects 
can be expected to lead to some events failing to satisfy 
(HI). In our simulations, the consistency requirement (|8]) 
rejects about 30% of the remaining events from both the 
signal and the Z'^ — >■ tt background sample. 

Example distributions for TT^H^gg^-bound ^^^^^^ g, 
number of other existing kinematical variables, described 
later) can be found at Monte Carlo truth level in Fig- 
ure [1] and after basic detector simulation in Figure [2] 
We recall that a perfect, hermetic, detector would guar- 
entee that TTjH^ggs-bound ^^j. ^.j^g signal and similarly 

^mggs-bound < fQj, ^i^g Z ^ TT background. When 
detector resolution effects are added (Figure [2|) we ob- 
serve a small tail of events creeping above the ideal bound 
due to the smearing of the missing transverse momentum. 
Nevertheless, the signal and background Tj^H^ggs-bound 
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FIG. 2: The distributions of the various discriminating variables, at detector level and after the selection cuts for the simulated 
120 GeV h ^ TT signal (solid) and Z -> rr background (dashed). The new variable m™'^'^""'""""'* 
transverse mass m^^'^ |(b)| and two other variables currently employed by the LHC collaborations, m 



IS shown, as is the 
and mTf^^" 



(d) All signal and background histograms are presented with unit normalization to allow comparison of shape. 
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FIG. 3: Optimal discovery potential (S/^/B) using each of 
the different discriminating variables. The best discrimina- 
tion is found using m!^^^^''~''°"°'* for those events where it 
exists and 771^"^"° for the remainder. 



imisation is reduced. This leads to a large fraction of the 
events lying close to the the upper kinematic end-point. 
The resulting distributions are then sharply peaked, and 
show good separation between the Higgs Boson signal 
and the dominant Standard Model background process. 

Figure O also shows unit-normalized distributions of 
three other mass-sensitive variables. The transverse mass 
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to rn^'J^gs-honnd 



shown in Figure 2(b) is defined in a similar manner 
but lacks the t mass-shell constraints 
Figure 2(c) shows the distribution of the so- 



called "tau-tau effective mass" , defined by 



m. 



EffectivcN ■' 



(Pf 



P2' 



P'. 



2ti 



where i?^ = {p^,pr^,0) is a massless four- vector con- 
structed from the missing transverse momentum. This 
variable has been used as a discriminant for fully leptonic 
(e^, /i^) tau events. In our simulations it generates broad 
distributions with rather poor separation between signal 



and background. Figure 2(d) shows the distributions for 



the invariant mass of the visible decay products 



distributions fall off rapidly above rrih and mz respec- 
tively. We note that because Tj^H^sss-bound j^jg^j^gg ^gg of 
the full set of kinematic constraints ([3H7]) , the space over 
the hypothesized momenta can be chosen during the min- 



(m™^)' = (PC 



P2niPi,+P2^), 



a discriminant that has been used in observation 
Z — )■ rr [1^ and searches for MSSM Higgs bosons P, | 



of 
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Again, the distributions are broader, and the separation 
between signal and background less pronounced, than for 

j^Higgs-bound 

We quantify the discrimination power of each of the 
four variables as a function of rrih as follows. For each 
variable / e {m™'='=*'™, my^^^^'^ ^^8^''-'"°""'^, to|™'=} 
we plot the distribution, after detector simulation and the 
cuts described above, of the signal and of the dominant 
Z — > TT background. We then determine the value of 
the additional cut that maximises S/^/B subject to 
S > 10, where S and B are the numbers of signal and 
background events, respectively, that would be expected 
with / > J™'" for 15 fb^^ of integrated luminosity. 

The optimal values of S/^/B obtained for the four dif- 
ferent projections are shown as a function of mh in Fig- 
ure O The line marked ^"^88""*^°""'^ is the result ob- 
tained if one simply discards the ^30% of events that 
fail to satisfy dH). and hence for which TO^^ss^-bound 
not defined. 

Both ml™" and m,™gg>5-bound perform better than 
the two currently-employed alternatives. We might ex- 
pect 777,H^gss-bound^ which has a sharper distribution than 
m!P'"°, to have a larger significance. However, because it 
is only possible to define TTjH^sss-bound -j-^j. gyents that sat- 
isfy ([8]), its stand-alone significance is reduced. What one 
really wants is to combine the desirable features of the 
two related mass-bound variables — the sharper distribu- 
tion of 77^H|ggs-bound^ ^^iQ guarcntccd existance, even 
after smearing, of to|'"°. One way of doing so is by form- 
ing a 'best effort' variable, defined to be ttjH^sss- bound ■£ 
it exists, and Tn^'^^° otherwise. The resulting distribution 
is found to have the largest value of S/ \fB for all values 
of nih simulated. 

The peformance of the Bound-or-True combined vari- 
able has been compared to yet another alternative, one 
obtained from an event- by- event maximization of a likeli- 
hood over all allowable neutrino momenta |2C|. The dis- 
tribution obtained using the method of Ref. [20| results 
in an almost equally high discovery potential, but unlike 
our proposal it must be tuned to the particular kinematic 
cuts empolyed, and it is a factor of >1000 more compu- 
tationally expensive to compute. 

Conclusions 

We have advocated use of TTjH^sgs-bound -j.^ separate 
h ^ TT events from their most significant irreducible 
Standard Model background Z — >■ tt. We observe that, 
for events in which a solution to the full set of kinematic 
constraints exists, the superior discriminatory power of 
^Higgs-bound a^j-jggg j^y construction: it is as the maximal 
lower bound on an important scale (in this case the in- 
variant mass of the tau pair) , which is smaller than mz 
for the background, but is often ~ mt for the signal. 

For events with no TOH^ggs-bound solutions, the usual 
transverse mass m^™'^ can be substituted in the same 



role, so that one makes good use of all events. Choosing 
j^Higgs-bound -^^j^gjj available and m^^° otherwise offers 
superior discriminatory power than either variable alone, 
and in all cases a better discovery potential than the 
variables currently being used in LHC searches. 

Although we have focused our attention here on the 
decay h — > tt, it is worth remembering that the same 
variable can be used on any topology with a similar kine- 
matic structure having a resonance decaying to a pair 
of intermediate on-shell particles under the substitution 

y ^intermediate- 
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